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Phase Structure of Lattice QCD with Wilson fermion at Finite 
Temperature 

S. Aoki^ 

'^Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 

We review our attempt at understanding the phase structure of lattice QCD with the Wilson fermion formu- 
lation at finite temperature, based on a spontaneous breakdown of parity-flavor symmetry. Numerical results 
demonstrating explicitly the spontaneous breakdown of parity-flavor symmetry beyond the critical line at both 
zero and finite temperatures are reported. The phase structure and order of the chiral transition are reported for 
Nf = 2,3, and 4 fiavors, and the approach towards the continuum limit is discussed. 
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Figure 1. A phase structure at finite temperature. 

1. Introduction 

Studying the finite-temperature QCD phase 
transition via the lattice regularization with the 
Wilson fermion formulation is known to be more 
difficult than with the staggered fermion formu- 
lation. This is because of absence of an order 
parameter due to explicit chiral symmetry break- 
ing introduced by the Wilson term. A naive ex- 
pectation for the phase structure in a /3 = 6/5^ 
vs. K = l/(2mga-|-8) plane at finite temperature 
with fixed temporal lattice size which was the 
basis of early studies of the phase transition is 
sketched in Fig. |l|. At zero temperature the pion 



mass vanishes along a line Kd/S), the critical line, 
which runs from Kc{P = 0) = 1/4 in the strong 
coupling limit to Kc{(3 — 00) — 1/8 in the weak- 
coupling limit. At finite temperature, which is 
realized at fixed Nt <^ L where L is a spatial 
lattice size, another critical line Kt{/3), the fi- 
nite temperature phase transition line, should ap- 
pear starting from K — 0, which corresponds to 
the deconfinement transition of pure gauge the- 
ory. The critical line, though rigorously speaking 
some part of the line may merely show cross-over 
phenomena, is expected to extend into the non- 
zero K region from below and hopefully to merge 
with the critical line Kc at some /3. The crossing 
point, if it exists, may be identified as the point 
of chiral phase transition, since at this point the 
finite temperature phase transition occurs in the 
presence of a massless pion, or equivalently mass- 
less quarks. 

In Fig. |l| the region on the left of Kt and be- 
low Kc is in the low temperature phase, while 
the region on the right of Kt is in the high tem- 
perature phase. The critical line Kc in the high 
temperature phase is no more critical, since the 
pion mass does not vanish in the high temper- 
ature phase, as naively expected and confirmed 
by numerical simulations. Therefore in Fig. |l| a 
dotted line is assigned to the part of the Kc line 
included in a high temperature phase, where the 
pion mass does not vanish. 

Although the phase structure above seems 
plausible at first sight, several questions arise. 
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First of all, what is the meaning of Kc at finite 
temperature ? We define Kc as a line of vanish- 
ing pion mass at zero temperature, but according 
to the phase structure above the pion does not 
become massless on Kc in the high temperature 
phase. This suggests that a location of the line 
Kc may depend on temperature if it is defined 
as a line of vanishing pion mass at correspond- 
ing temperature. A related question is the na- 
ture of the region above Kc at small /3, where a 
question mark is assigned in Fig. ^ This region 
seems to belong to the low temperature phase 
since the pion can become arbitrarily light in the 
region near Kc- If this is the case, the line Kt 
should continue to exist beyond the chiral tran- 
sition point in larger K region, so that the low 
temperature phase above the Kc can be separated 
from the high temperature phase. This should 
be numerically confirmed. A more fundamental 
question is how one can understand the existence 
of a massless pion at zero temperature without 
recourse to chiral symmetry. 

In this talk we try to answer these questions, 
based on interpretation of the massless pion in the 
Wilson fermion formulation, using both analytic 
methods and numerical simulations. 

2. Phase structure at zero temperature 

2.1. Parity-flavor breaking scenario 

In this subsection we briefly review our under- 
standing of massless pions appearing at Kc in the 
Wilson fermion formulation at zero temperature 
based on refs.|^-||. 

Let us consider 2- flavor QCD. Since pions be- 
come massless dX K = Kc we expect the following 
behavior of the 2 point function for the corre- 
sponding operators tt"" — tpij^r'^ip {a — 1, 2, 3): 

hm l^(7r'^(0,0K(x,0) 

X 

r ZiexpH/^], K<Kc 

= { Z2t-", K^Kc 

[ (7r«)2 + Z3CxpH/e], K>Kc 

where the correlation length ^ is related to the 
mass of the vr meson as ^ = l/{mT^a), a is some 
exponent, and V = is the space volume. 



Since an existence of massless tt" meson implies 
a divergent correlation length at Kc, the opera- 
tor corresponding to the tt" meson is expected to 
develop a long-range order for K > Kc so that 
{tt"') 7^ 0. This non-zero condensation is possible 
for all a in general, but using flavor symmetry 
we can make a change of variables such that only 
one tt" develops a non-zero vacuum expectation 
value. Hereafter, without loss of generality, we 
assume that (tt^) ^ and (tt^^^) = 0. 

This condensation spontaneously breaks both 
parity and flavor symmetry, and a second order 
phase transition from the parity-flavor symmet- 
ric phase to a parity-flavor broken one occurs 
at Kc- Therefore a massless (neutral) tt me- 
son corresponding to the tt^ operator can be in- 
terpretated as a massless mode associated with 
this second order phase transition, not a Nambu- 
Goldstone boson. This neutral tt meson becomes 
massless only at the transition point K = Kc, 
while the charged tt mesons corresponding to the 
linear combinations of tt^ and tt^ operators stay 
massless in the broken phase since they are the 
Nambu-Goldstone bosons associated with the fla- 
vor symmetry breaking. 

According to the above interpretation we ex- 
pect near Kc, {7x1^,0)"^ ~ [Kc - K^" oc (M - 
Mc)'^'^ with a critical exponent i', where we de- 
fine M = 1/{2K) for later uses. Since low en- 
ergy properties of tt mesons may be described by 
an effective 4-dimensional scalar field theory, we 
can expect the phase transition to be mean-field 
like up to logarithmic corrections and therefore 
1/ = 1/2, which reproduces the well-known PCAC 
relation {m-^a)^ oc rrigO where quark mass may be 
defined as m^a ~ M — Mc- 

On the other hand, if we define the current 
quark mass via the axial Ward identities]^: 

2^wi ^ Aa;(0|Vii757oV'(a^)k) 

this automatically satisfies the PCAC re- 

lation cx at small mj. However, since 

this rn^^ is not a tunable parameter, an existence 
of the chiral limit where 711-^ becomes zero is not 
ensured. Our scenario explains why such a tun- 
ning is possible. Note that (x [M — Mc) if 
V = 1/2. 
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Figure 2. (a) {ipi^BT^ip) as a function oi M — 
1/2K at /? = on a 4^ lattice with H = 0.05, 0.1, 
0.2, 0.3 (solid symbols) and on a 6'' lattice with 
H = 0.05 (open circles). The lines are analytical 
predictions obtained at /? = in the large Nc limit 
with the corresponding H. (b) same for (ipip) . 



This new scenario for the existence of the mass- 
less pion in the Wilson fermion formulation has 
been confirmed analytically in the strong coupling 
and expansions. See refs.[^|,^ for details. 

2.2. Numerical evidence for the parity- 
flavor breaking phase 

In this subsection numerical evidence will be 
given for the existence of the parity-flavor break- 
ing phase in full QCD with Wilson fermions. 

In order to examine if tt^ = iipjsT^ip can 
develop a non-zero expectation value, we have 
to add a source term iH'tp'^^T^ip to the Wilson 
fermion action. The reason for this is as follows. 



Since Re[z^75r'^^] changes its sign according to 
the sign of H for an arbitrary configuration, it is 
easy to show, at = on a finite lattice, that 

Re[ii}^5T^i}] = 

without averaging over configurations. Therefore 
one can not see any signal of a non-zero value of 
{i-tfj^^T^tjj) from simulations with H = 0. This 
is in contrast to the case of the magnetization 
of the Ising model, where it takes non-zero val- 
ues on a given configuration for zero magnetic 
field and may change its sign from configuration 
to configuration. To obtain the correct value of 
{iip^^T^ip) in Wilson fermion formulation, one 
has to take limits in the following order: 

lim lim (i^j^T^^). 

This means that one should make an extrapola- 
tion to -ff = from a series of values at non-zero 
H's on finite but sufficiently large lattices. 

Fortunately, even in the presence of the source 
term, the Wilson fermion determinant is still real 
positive, so that the hybrid Monte Carlo method 
can be utilized. 

In order to establish the parity-flavor breaking 
phase we have carried out a full QCD simulation 
with Nf = 2 at /3 = 0, where analytic expressions 
for the order parameters are available for non-zero 
external field H within 1/Nc expansion 

In Fig. ^, (ipij^T'^ip) and {"ipip) are plotted as a 
function of M = 1/2K at /3 = on a 4"* lattice 
with H = 0.05,0.1,0.2,0.3 and on a 6'' lattice 
with H = 0.05, together with the analytic pre- 
dictions with the corresponding H in the strong 
coupling and large Nc limits. Since the agreement 
between the numerical data and the analytic pre- 
dictions is very good, it can be established at least 
in the strong coupling limit that the parity and 
flavor symmetry breaking phase transition really 
occurs at M = Mc = 2 {K = Kc = 1/4) and 
{tpij^T^ip) is non-zero in the iJ — > limit at 
M < Mc {K > Kc). We have also shown that 
finite size effects on the order parameter is very 
small by comparing data on a 4^ lattice with those 
on a 6* lattice or with the analytic predictions 
which correspond to the values in the infinite vol- 
ume limit. 
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Figure 3. Susceptibility X55 as a function of M — 
1/2K at /? = on a 4^ lattice (solid symbols) and 
on a 6^ lattice (open circles). 



As a non-trivial check we have calculated sus- 
ceptibilities X55: XoOj and xso; which are defined 

by 

77 ( i'nuT^^ix) ) 



X55 



Xoo 



X5Q 



\ X y ^ 

We have plotted X55 in Fig. | and xoo, Xso in 
Fig. ^, together with the analytic predictions. 
The agreement between the numerical data and 
the analytic predictions is reasonably good, ex- 
cept for deviations observed deep in the bro- 
ken phase, where our analytic formula becomes 
inaccurate in the presence of massless Nambu- 
Goldstone charged pious. These results for the 
susceptibilities also support our conclusion about 
the existence of the parity-flavor breaking phase. 

2.3. Phase structure at small p 

In the previous subsection we have numerically 
established that the parity-flavor breaking phase 
indeed exists in the strong coupling limit. This 
phase should continue to exist in the strong cou- 
pling region, since a strong-coupling expansion, 
which support an existence of the phase |3|, is con- 
vergent around /3 = 0. However an understand- 
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Figure 4. (a) xoo- (b) xso- 



ing of the phase structure in the weak coupling 
region is particularly important, since the con- 
tinuum limit can be defined only in this region. 
Therefore, in this subsection, we consider how the 
phase structure changes toward the weak coupling 
region. 

In Fig. ^ we propose the expected phase struc- 
ture of lattice QCD with 2 flavor Wilson fermions 
in the jS vs. AI = 1/2K plane, based on the rea- 
soning mentioned later. At (3 < (3^ there are only 
two critical lines, while there are ten critical lines 
at /3 > Puj, and 5 points where two lines meet 
at /3 = oo. At M = 4 on can define the usual 
continuum limit where only a physical mode of 
the sixteen zero modes of the lattice fermion be- 
comes light, while the continuum limits at other 4 
points contains light doubler modes. For example 
4 doubler modes, which has one of the momentum 
component near ir/a, appear as light fermions at 
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Figure 5. The expected phase structure of lattice 
QCD with Wilson fermion in /? - M plane. 



M = 2. Note that the phase diagram is shown to 
be symmetric under the M — > — M transforma- 
tion. 

There are several evidence for the existence of 
multiple critical lines in the weak coupling region. 
First of all, a phase structure similar to Fig. || has 
been found for the Gross- Neveu model in 2 dimen- 
sions in the large N limit, regulated on the lattice 
using Wilson fermions . The existence of a 
multiple structure of the critical lines has been 
also suggested from studies of eigen- value spectra 
of the Wilson fermion matrix on quenched config- 
urations!^,^ and on full QCD configurations]^. 

To establish the multiple structure of the crit- 
ical lines we performed new simulation using 
quenched configurations without adding the ex- 
ternal field iH^l^jii^p. We have calculated the pion 
mass at /? = 6.0 on a 10^ x 20 and a 16^ x 20 lat- 
tice to find the location of multiple K^s. The 
pion propagator is averaged over 10-20 configu- 
rations. 

We plot our results for the pion mass squared in 
Fig. ^ as a function of M ~ 1/2K. A reasonable 
agreement of results for lO'^ and 16'^ spatial sizes 
shows that finite spatial size effects are not severe 
in our data. We observe that there exist four more 
critical values of vanishing pion mass in addition 
to the conventional one located at the right most 
of the figure. The existence of 5 critical points for 
M > at fixed /3 supports the phase structure 
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Figure 6. Pion mass squared at zero temper- 
ature as a function of M = 1/2K at /3 = 6.0. 
Open circles are for a 10"^ x 20 lattice, and filled 
circles for a 16^ x 20 lattice. Lines show linear 
fits of in M. 



proposed in Fig. ^. 

3. Phase structure at finite temperature 

Since we now understand the phase structure of 
lattice QCD with the Wilson fermion formulation 
at zero temperature, we can extend our investi- 
gation to the case at finite temperature. 

3.1. Expected piiase structure at fixed Nt 

In this subsection, based on the finding at zero 
temperature, we propose the phase structure for 
a fixed temporal extension Nt corresponding to 
finite temperature. 

We first consider how the multiple structure 
of critical lines is modified for a finite Nt- An 
important hint for the answer to this question 
again come from the Gross-Neveu model in 2 di- 
mensions. In Fig. we plot the critical line cor- 
responding to TOtt ~ on the {g, m) plane for 
Nt = 00,16,8,4 and 2 where g and m are the 
bare coupling and quark mass, respectively, the 
latter related to the hopping parameter through 
K = l/(2m -I- 4). The crucial feature revealed in 
Fig. ^ is that the three cusps of the critical line 
retract from the weak coupling limit g = for 
finite temporal lattice sizes, forming a single con- 



6 




1.5 1.0 0.5 

g 



Figure 7. Critical line for the lattice Gross-Neveu 
model on [g, m) plane. Temporal lattice size 
equals Nt — oo, 16, 8, 4, and 2 from outside to 
inside. 
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tinuous line which shifts toward strong coupling 
as Nt decreases. Thus, for a finite Nt, the criti- 
cal line is absent for sufficiently weak coupling. 

The close resemblance of the Gross-Neveu 
model and QCD such as the asymptotic freedom 
and the spontaneous chiral symmetry breaking 
suggests that a similar behavior of the critical 
line will hold for QCD, except that the number of 
cusps should increase to five because of the differ- 
ence in the number of dimensions, as seen in the 
previous section. For the finite temperature case, 
an apparent disappearance of the critical line to- 
ward weak coupling has been noticed in previ- 
ous studies jl^jl^. This means that Kc at finite 
temperature differs from Kc at zero temperature, 
if we define the Kc as a line of vanishing pion 
(screening) mass at the corresponding tempera- 
ture. 

In addition to the deformation of the critical 
lines at finite temperature, the thermal transition 
line Kt{P), which is absent in the Gross-Neveu 
model, should appear in QCD. As discussed in the 
introduction, in order to define the unique point 
of chiral phase transition at fixed Nt, the thermal 
line Kt{P) has to cross the critical line KdP). 
However the line Kt can not go into the parity- 
flavor breaking phase, since the massless pions 
exist in the phase so that it belongs to the low 
temperature phase. Therefore the line Kt never 



cross the line Kc defined at fixed Nt- In other 
words the region close to the critical line has to be 
in the low temperature phase even after it turns 
back toward strong coupling. This means that the 
thermal line should run past the turning point 
of the critical lineQ and continue toward larger 
value of K, as sketched in Fig. ^. Although Kt 
seems to cross Kc{T — 0) as seen in Fig. ^, no 
special criticality associated with massless quarks 
occurs at the crossing point. Therefore, if the 
above phase structure is true, we can not assign 
an unique point to define the point of chiral phase 
transition. 

3.2. Results for A^/ = 2 at At=4 

In this and the following subsections we review 
results supporting the phase structure proposed 
in the previous subsection. Here we report results 
for the 2-flavor case obtained for the temporal size 
Nt = 4 with the standard one plaquette action 
for gauge fields. 

We first show evidence that the cusps of the 
critical line retract from the weak coupling limit 
for a finite Nt as speculated in the previous sub- 
section. For this purpose we have carried out hy- 
brid Monte Carlo simulations of Nf — 2 full QCD 
using the Wilson quark action with = on an 

^Logically it may be possible that the hne Kt touches 
KciT 7^ 0) without crossing it. 
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Figure 9. tt and p screening masses and quark 
mass for N j = 2 full QCD with Wilson quark 
action obtained on an 8'^ x 4 lattice periodically 
doubled in one of spatial directions. Lines are 
linear fits to three data points. 



8'^ X 4 lattice] 12|. In Fig. ||(a) we present our re- 
sults for the TT and p screening masses, measured 
on the periodically doubled lattice in the spatial 
directions 8^ x 4 -> (8 x 2) x 8^ x 4, at /3 = 3.5 
as a function of \/K. Also plotted is the quark 
mass defined via the axial Ward identity, eq. (jl|) . 

As usual the pion mass squared and the quark 
mass linearly vanish almost simultaneously at the 
conventional critical point Kc ~ 0.2267, while 
another critical point at ~ 0.2454 can be 
clearly observed from the behavior of the pion 
mass squared and the quark mass in the larger 
K region. The existence of two critical points 
close to each other is quite analogous to Fig. ^ in 
quenched QCD at zero temperature. On the 

■^Although we have not carried out further simulations at 
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Figure 10. Parity- flavor order parameter as a 
function of 1/ A' for various values of external field 
H at (3 = 3.5 on an 8^ x 4 lattice for Nf = 2. 



other hand Fig. ^(b) shows how the behavior of 
these masses changes at /3 = 4.0. We observe that 
the gap between the two critical values has either 
become extremely narrow or disappeared. These 
results at two (3 values are consistent with the ex- 
pected deformation of the critical line at fixed Nt 
as shown in Fig. ^. 

Now one may ask whether the parity-flavor 
symmetry is indeed spontaneously broken in the 
region between two K^s a.t (3 = 3.5. In order to 
confirm this we have carried out new simulations 
with non-zero ij]!^. In fig. |l^ we plot the parity- 
flavor breaking order parameter < 'ipij5T3ip > as 
a function of l/K for H = 0.2,0.1,0.05,0.02 and 
0.01. Vertical lines mark the position of the crit- 
ical points at /3 = 3.5. We clearly observe that 
the order parameter between the two critical val- 
ues Kc and AT^ tends to a non-vanishing value as 
the external field H is reduced, while it decreases 
toward zero outside, supporting the spontaneous 
breakdown of parity-flavor symmetry between the 
two critical points. 

We next examine if the thermal line Kt run 
past the turning point of the critical line around 
P = 4.0 and continue towards larger values of 
AT, so that the region close to the critical line is 
in the low-temperature phase even after it turns 
back toward strong coupling. For that purpose we 



larger K to identify 3 more critical points, there is no 
reason to doubt their existence. 
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Figure 11. Real part of Polyakov line(i?eil), 
gluon (sg) and quark (sq) entropy density in lat- 
tice units for Nf = 2 full QCD with Wilson quark 
action obtained on an 8^ x 4 lattice. 



have measured thermodynamic observables, the 
real part of the Polyakov line and quark and gluon 
entropy densities in lattice units, and plot them 
as a function of l/if at /? = 3.5 in Fig. 



ll[a). 



The two vertical lines again show the position 
of the critical points estimated from the pion 
screening mass. Large values of these observables 
suggest that the system is in the high tempera- 
ture phase for smaller values of 1/K. However 
they decrease as 1/K increases toward the crit- 
ical point at 1/Kc « 4.075, becoming roughly 
similar in magnitude to those on the other side of 
the conventional critical point at 1/Kc ~ 4.411. 
Since the conventional critical point is in the low- 
temperature phase, a similarity of these thermo- 
dynamic quantities in magnitude support our ex- 
pectation that the other critical point is also in 
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Figure 12. Phase structure oi Nf — 2 QCD with 
Wilson quark action. Estimates of the critical line 
for Nt = 4 (solid line) and the thermal transition 
line (dotted line) are shown. 



the low-temperature phase. 

Behaviors of thermodynamic observables at 
f3 = 4.0 are shown in Fig. |ll|(b), together with the 
vertical line which marks the point where the two 
linear extrapolations of the pion screening mass 
squared in Fig. ^(b) cross each other. Since the 
thermodynamic quantities around the line stay 
small, the region around the "crossing point" is 
likely to be also in the low-temperature phase. On 
the other hand the increase of the three quantities 
across the line suggests that the thermal line runs 
very close to the critical line at /3 = 4.0. Within 
the resolution of our numerical simulation, how- 
ever, we can not exclude the possibility that the 
thermal transition line touches the critical line 
around P — A.O and run away again for larger K. 

Overall results in this subsection are consistent 
with our proposed phase diagram in Fig. ^. For 
definiteness, in Fig. |l^, we draw the phase dia- 
gram of two flavor QCD for Nt = 4, mainly based 
on our finding. We estimate the location of the 
critical line Kc (solid line) from the points of van- 
ishing pion screening mass (solid circles) as usual, 
supplementing the points where the inverse of the 
number of conjugate gradient iterations in the 
runs is linearly extrapolated to zcro(opcn circles). 
Here we assume that the turning point of the crit- 
ical line is located a,t /3 — 4.0. The dotted line rep- 
resents the thermal transition line Kt estimated 
from previous works (open circles) PJl^,|ll|. Our 
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Figure 13. tt screening mass squared and quark 
mass as a function of 1/K for Nf = 2 full QCD 
with RG improved gauge action and Wilson quark 
action at (a) (3=1.2 and (b) /3=1.35, obtained on 
an S'^ X 4 lattice periodically doubled in one of 
spatial directions. 



drawing of the dotted line continues above the 
upper part of the cusp of Kc, though its precise 
location can not be determined from our numer- 
ical data. 

3.3. RG improved gauge action for Nf—2 
3.3.1. Evidence of the cusp structure of the 
critical Hne 

From our data in the previous subsection the 
thermal transition around (3 = 3.5 ~ 4.0 seems 
to be a smooth cross-over at Nt = 4, and this 
agrees with the prediction that the chiral phase 
transition is of second order for Nf = 2|14| if the 



effect of the chiral anomaly remains at T « Tc. 
However it has been reported that the thermal 
transition becomes strongly first order around 
P = 5.0|l^, where the thermal transition line 
comes very close to the zero temperature criti- 
cal line. According to the phase diagram in the 
previous subsection, however, the critical line Kc 
defined for fixed Nt does not exist in the region 
where the first order transition is observed. In 
fact the pion screening mass never becomes small 
in the region of the first order transition. There- 
fore this first order transition has nothing to do 
with the chiral phase transition. Instead it is 
likely that it is a lattice artifact, and this view is 
indeed supported ||l5| by the fact that the strong 
first order transition disappears if one uses the 
renormalization group (RG) improved action for 
gauge fields p6[. 

One may wonder whether our findings such 
as the parity-fiavor breaking phase, the multiple 
structure of the critical lines, the formation of the 
cusp of the critical lines at finite temperature, and 
the relation of the thermal transition line to the 
critical line, are universal, and not a lattice arti- 
fact like the first order transition at /3 w 5.0 . 

To partly answer this question, we have made a 
new run using the RG improved gauge action with 
2 fiavors of Wilson quark action[|l^ . In Fig.|l^ the 
pion screening mass squared and the quark mass 
are shown as a function of l/-fir at (a) /3 = 1.2 
and (b) (3 = 1.35 on an 8'^ x 4 lattice, periodi- 
cally doubled in one of spatial directions for mass 
measurements. These results show that there are 
two critical points K^ at /3 = 1.2, whose gap then 
becomes very narrower or disappears at /3 = 1.35. 

Therefore at least the multiple structure of the 
critical lines and the formation of the cusp of 
the critical line at finite temperature are common 
properties in both the standard gauge action and 
the RG improved one, and they may well be not 
a lattice artifact. 

3.3.2. ScaUng analysis 

Assuming the chiral phase transition for 2 fla- 
vor QCD is of second order, the magnetization 
M is expected to be described by a single scaling 
function, 
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Figure 14. Best fits to the scaling function with 
(a) 0(4) and (b) MF exponents. The plots con- 
tain all the data within the range < 2mqa < 0.8 
and p < 2.0. Solid curves are scaling functions 
obtained in an 0(4) spin modeljl^] and in a MF 
calculation, respectively. 



where h is the external magnetic field and t — 
[T - Tc{h = 0)]/Tc{h = 0) is the reduced tem- 
perature. For 3 dimensional 0(4) models, the 
critical exponents are given by l/(/3<5) = 0.537(7) 
and 1/(5 = 0.2061(9) @. 

Due to the explicit chiral symmetry breaking 
of the Wilson quark, eq. @) contains the finite 
lattice spacing errors, which may be reduced by 
the improved action. Indeed it has been reported 
that eq. (^ is numerically satisfied for the RG 
improved action |l^ as follows. 

Usually M is identified as {"ipip)- For Wilson 
quarks, a proper subtraction is necessary to ob- 
tain the correct continuum limit. The properly 
subtracted (ipip) defined via the axial Ward iden- 



tity (| as 

{^4>)suh = 2TO,a(27^)2 ^(^(x)7r(0)) (3) 

X 

is proposed to be used for M, together with 
h ~ 2niq and t — [3 — Pet, where (3ct is a free 
parameter Assuming that the scaling func- 
tion f{x) is universal, and thus using the one ob- 
tained for an 0(4) model |jl^, one can fit the 
data by adjusting (3ct and the scales for t and 
h, with the exponents fixed to the 0(4) values. 
This works well with (3ct = 1-35, as shown in 
Fig. 0(a). Since it is suggested that the QCD 
chiral transition might be described by the mean- 
field (MF) criticalitypol, one may also test if the 
MF scaling function together with the MF expo- 
nents, l/{/36) = 2/3 and 1/6 — 1/3, is consistent 
with the data. Fig. ^(b) shows, however, that 
it gives no reasonable fit||l^. The success of this 
scaling test with 0(4) exponents suggest not only 
that the chiral transition in the continuum limit 
is described by the 0(4) criticality but also that 
the chiral transition is of second order for 2 fla- 
vors. It also shows that the RG improved gauge 
action indeed reduces the violation of the scaling 
relation, eq. (|^). 

3.4. Results for 3 and more flavors at Nt = 

4 

The thermal chiral transition is a smooth cross- 
over for 2 flavors of Wilson quark action, and 
hopefully it will show the criticality of the second 
order phase transition in the continuum limit. On 
the other hand, the chiral transition of QCD for 
Nf > 3 is predicted to be of flrst order in the 
continuum p^. Therefore, the thermal transition 
with 3 or more light quarks should also be of first 
order at least in the scaling region, even though 
the chiral symmetry is explicitly broken by the 
Wilson term. 

Actually a strong first-order transition has 
been observed for Nf = 3 at /3 « 4.0 ~ 4.7 on 
an 8^ X 10 X 4 latticejH]. Assuming the behav- 
ior of the critical lines at fixed Nt is similar for 
Nf = 2 and 3, their data also tell us that the 
tip of the cusp of the critical line, if it exists, is 
located around /? « 3.0. Like the Nf = 2 case 
the signal for the first order transition is only ob- 
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Figure 15. tt screening masses and quark mass as 
a function of l/K for A^/ = 3 full QCD with the 
Wilson quark action obtained on an 8^ x 10 x 4 
lattice periodically doubled in a spatial direction 
of length 10. 
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Figure 16. Polyakov loop as a function oi 1/K 
along the lines of K = Kcif3), K = K^{fi) - 0.005 
and K = K^{I3) - 0.01, for iV/ = 3 full QCD on 
an 8^ X 10 X 4 lattice. 



served far away from this tip of the cusp. There- 
fore, in order to claim that the thermal transition 
with three light Wilson quarks is of first order 
at Nt = 4, one should find the first order signal 
near the tip of the cusp around (5 w 3.0. A more 
fundamental question is whether the phase struc- 
ture we find for Nf = 2, such as the multiple 
structure of the critical lines and the formation 
of the cusp at fixed Nt holds also for Nf ~ 3 
or more. Moreover if the thermal transition near 
the tip of the cusp is of first order as expected, we 
may determine the precise location of Kt, so that 
the question whether the thermal line Kt and the 
cusp of the critical line Kc touch each other or 
not may be clarified, whereas it is difhcult to an- 
swer it for a smooth cross-over like in the Nf=2 
case. 

In order to answer these questions we have 
made a new run for Nf = 3 QCD using the Hy- 
brid R algorithm [itI]. In Fig. ^ the pion mass 
squared and the quark mass are again plotted as 
a function of 1/K at /3 = 2.5 on an 8^ x 10 x 4 
lattice. The linear extrapolation of these quanti- 
ties confirms the existence of two Kc's. Like the 
Nf = 2 case the multiple structure of the criti- 
cal lines presents for Nf = 3. Although we have 
not made similar measurement at (3 — 3.0, we 
expect, from the data of previous studyElt], that 



the critical line turns back toward the strong cou- 
pling limit, forming a cusp around (3 — 3.0. 

As far as the finite temperature chiral transi- 
tion is concerned, we observe that the strong first- 
order signal previously found around /3 = 4.7 be- 
comes weaker toward the stronger coupling, as is 
in the case of Nf = 2. To find a signal of the first- 
order transition near the tip of the cusp, we have 
made two sequences of runs varying K and f3 in 
accordance with relation that K — KdP) — 0.005 
OT K = Kc{f3) - 0.010, where K^ifi) is the crit- 
ical line defined at T = 0, not at finite temper- 
ature. Two state signals very close to the tip 
of the cusp, which is expected to exist around 
K = 0.235(: 1/K = 4.255 ) and /3 = 3.0, are 
observed. In Fig. |l^, we plotted the behavior of 
the Polyakov loop as a function of 1 /K along the 
lines oiK = Kc{f3) - 0.005 or K = K^{(3) - 0.01, 
together with the data previously obtained on the 
line of K c{f3)\2:l\{^ where the left most point may 
not be thermalized) . The signals of the first-order 
transition are seen. Therefore the thermal tran- 
sition around the tip of the cusp at Nt = 4 may 
be of first order for 3 flavors, as expected from 
the continuum. It is noted, however, that this 
data is taken for fixed step size 5t = 0.01 of the 
Hybrid R algorithm. Therefore more runs with 
a smaller 5t is need to definitely claim the exis- 
tence of the first order phase transition near the 
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for Nt = 4. The previous result for Nf = 2 is 
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signals are found, while open rectangles represent 
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Figure 18. The expected phase structure for dif- 
ferent Nt^s. 



tip of the cusp. 

Now we increase the number of flavors to 4. 
One of the advantage for Nf=A over the Nf = 3 
case is that the HMC method can be directly ap- 
plied to it, and thus no extrapolation of St is nec- 
essary. Our HMC runs for Nf = 4 with H = are 
carried out on an 8'^ x 4 lattice, and hadron masses 
are calculated on a periodically doubled lattice in 
one of spatial directions . The location of the 
critical line is estimated from the screening pion 
mass squared and the position and order of the 
thermal transition is examined through the be- 
havior of physical quantities. 

In Fig. ^ we show the phase diagram for 
Nf = 4: together with that for TV/ = 2 in Sec. |3^. 



We find again a cusp structure similar to that of 
Nf = 2 except for a shift oi 6/3 « 1.2 towards 
stronger coupling, which we qualitatively expect 
from a larger magnitude of sea quark effects for 
Nf — 4. We also find strong first-order signals 
across the thermal line away from the tip of the 
cusp as marked by solid squares. Similar to the 



cases of N 



f 



2 and 3, the transition becomes 



weaker towards the tip of the cusp, apparently 
turning into a smooth crossover at /3 = 2.5 — 2.0. 
The location of the crossover is indicated by open 
rectangles in Fig. nn 



One surprising thing is that no signal of the 
first order transition can be found so far near the 
tip of the cusp, contrary to our expectation. New 
runs are also made with an increased spatial size 
of 12"^ X 4. Results do not show any deviation from 
those on an 8'^ x 4 lattice. Therefore it is unlikely 
that finite-size effects have rounded a first-order 
discontinuity into a smooth crossover on an 8'^ x 4 
lattice. To find the signal of the expected first or- 
der phase transition near the tip of cusp, as in 
the case oi Nf = 3, simulations on the lines of 
K = Kc{l3) - 0.005/0.01 should be performed. 
We reserve, however, this study for future inves- 
tigations. 

4. Approach toward the continuum limit 

It is usually though that the continuum physics 
of finite temperature QCD can be extracted from 
the lattice simulations by the following procedure. 
For fixed iV^ the chiral transition point may be 
defined at the crossing point of Kc and Kt, and 
physical observables are measured at that point. 
One should repeat this for increased Nt, which 
corresponds to the decreased lattice spacing as 
a 1/Nt, to extract the continuum values of the 
observables. 

Based on our understanding of finite temper- 
ature phase structures with Wilson fcrmion for- 
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Figiire 19. Location of the critical line for the 
lattice sizes 8^ x 4 (dotted lines) and 8^ x 8 (solid 
lines) for Nf = 2 and 4. 



mulation, however, it now becomes clear that the 
thermal transition point at fixed Nt in the pres- 
ence of light quarks can be defined anywhere in 
the region around the tip of the cusp. There- 
fore, in principle, no unique definition of the chi- 
ral transition point exists at fixed Nt- In prac- 
tice, it is better to choose one definite prescrip- 
tion for the chiral transition point. For example, 
as shown in Fig. |l^, the lower crossing point of 
KriNt) and Kc{T = 0) may be regarded as the 
chiral transition point, on which observables such 
as Tel nip, etc have to be measured. To approach 
the continuum limit and to investigate the scal- 
ing behavior, one has to increase Nt, repeat same 
measurements on the chiral transition point cor- 
responding to increased Nt- 

One may choose freely another definition for 
the chiral limit, as the difference between two 
choices is 0{a) and should disappear in the con- 
tinuum limit. For example, the chiral limit may 
also be defined at the crossing point of KT{Nt) 
and the line Kq {T) where the quark mass rri^'^ of 
eq-dl) vanishes. Note that Kc{T = 0) = Kq{T ^ 
0) but Kc{T T^O) ^ Kq{T ^ 0) in general at fi- 
nite temperature. In particular Kq(T) seems to 
run from the low temperature phase into the high 
temperature phase, so that it can cross KT{Nt). 

In order to approach the continuum limit in 
this way, it is important to know how the cusp of 
the critical line moves when the temporal lattice 
size is increased. To get a hint for an answer 



to this question, we estimate the location of the 
critical lines on an 8'^ x 8 lattice for Nf — 2 and 
4. 

In Fig. |l^ we plot the critical line estimated 
from the pion mass squared. Although the cusp 
is shifted toward weak coupling, the magnitude 
of shift is quite small both for Nf = 2 and 4, so 
that the tip of cusp is still in the region of strong 
coupling for Nt = 8. The slow move of the tip has 
already been indicated in the previous studies for 
Nf — 2 IH], whose data can be interpreted that 
the tip of the cusp is located at /3 ~ 4.0 — 4.2 for 
A^f = 6 and at /3 = 4.5 - 5.0 even for Nt = 18. If 
this is indeed the case, a substantial increase in 
the temporal lattice size is needed for the tip of 
the cusp and the thermal transition line moving 
into the scaling region, which is usually though 
to start at /3 « 5.5 for Nf = 2. 

It should be noted that the slow move of the 
tip of the cusp has a significant impact also on 
spectrum calculations at zero temperature. Since 
the location of the tip of the cusp is controlled 
dominantly by the smallest extension of the lat- 
tice, the very large spatial size is necessary for an 
existence of the critical hopping parameter, thus 
an appearance of massless pions, in the scaling re- 
gion. Otherwise the finite lattice size effect may 
become very large near the chiral limit at larger /3, 
or to say differently, it could happen that hadron 
masses are measured at some /3 where the critical 
line is absent. 

5. Conclusions 

In this talk we have tried to present our un- 
derstanding on the phase structure of QCD with 
Wilson fermion both at zero and finite temper- 
atures, based on the spontaneous breakdown of 
parity-flavor symmetry. We have shown that so 
far all existing numerical data as well as analytic 
calculations are consistent with our understand- 
ing. 

In future study it is intersting to investigate 
the phase structure of QCD with the improved 
Wilson fermion action (: clover action). Since we 
believe that the parity-flavor breaking mechanism 
to obtain massless pions is quite general for the 
Wilson fermion type action who looses the chiral 
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symmetry, a similar phase structure is expected 
to exist. Furthermore it is important to see how 
fast the tip of cusp, if it exists, moves toward 
the scahng region as Nt increases. Hopefully the 
improvement of the fermion action as well as that 
for the gauge action may overcome the problem 
of the slow move in the previous section. 

Although a lot of things concerning the parity- 
flavor breaking mechanism still remain to be in- 
vestigated, we think that no conceptual difficulty 
in Wilson-type fermion actions at finite temper- 
ature exists. So now it is a good time to start 
serious study, like the Kogut-Susskind quark ac- 
tion, to get real physical results from the lattice 
QCD with Wilson-type fermion action, but it may 
be necessary to use gauge and fermion improve- 
ments @^|2). 
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